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Optimal design of PID controllers using the QFT method 
A. C. Zolotas 1 and G. D. Halikias 2 



Abstract: An optimisation algorithm is proposed for designing PID controllers, 
C\f- which minimises the asymptotic open- loop gain of a system, subject to appropriate robust- 

stability and performance QFT constraints. The algorithm is simple and can be used to 
automate the loop-shaping step of the QFT design procedure. The effectiveness of the 



method is illustrated with an example. 
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1 Introduction 



Many practical systems are characterised by high uncertainty which makes it difficult to 
O ■ maintain good stability-margins and performance properties for the closed-loop system. 



There are two general design methodologies for dealing with the effects of uncertainty: 
^. ■ (i) Adaptive control, in which the parameters of the plant (or some other appropriate 

structure) are identified on-line and the information obtained is then used to "tune" the 



controller, and (ii) Robust control, which typically involves a "worst-case" design approach 
for a family of plants (representing the uncertainty) using a single fixed controller. 

Quantitative Feedback Theory (QFT) is a robust-control method developed during the 
last two decades which deals with the effects of uncertainty systematically. It has been 
sucessfully applied to the design of both SISO and MIMO systems, while the theory has 
^ " also been extended to the nonlinear and the time-varying case. In comparison to other 

optimisation-based robust control methods, QFT offers a number of advantages. These 
include, (i) the ability to assess quantitatively the "cost of feedback" [3], [4], [5], (ii) the 
ability to take into account phase information in the design process (which is lost if, e.g. 
singular values are used as the design parameters) , and , (iii) the ability to provide "design 
transparency" , i.e. clear tradeoff criteria between controller complexity and feasibility of 
the design objectives. Note that (iii) implies in practice that QFT often results in simple 
controllers which are easy to implement. 

For the purposes of QFT, the feedback system is normally described by the two- 
degrees-of-freedom structure shown in Figure 1. In this diagram, G(s) is the uncertain 
plant, K(s) is the feedback controller and F(s) is the pre-filter. The objective is to design 
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K(s) and F(s) so that the output signal Y(s) tracks accurately the reference signal R(s) 
and rejects the disturbance D(s), despite the presence of uncertainty in G(s). 

D(s) 




A - 



Figure 1: Two degree of freedom feedback system 

The uncertain plant G(s) is assumed to belong to a set, G(s) G {G(p, s) : p G V}, 
where p is the vector of uncertain parameters; these are assumed to be fixed but unknown, 
and to take values in V. The uncertain plant is first translated in the frequency domain 
(using a discrete grid of frequencies u>i, 0J2, ojn, typically chosen to cover adequately the 
system's bandwidth), resulting in N "uncertainty templates" Qi = {G(p,jUi) : p G V}, 
i = 1,2,..., TV. The tracking specifications are given in the form of upper (B u (ui)) and 
lower (Bi(ui)) bounds in the frequency domain, usually based on simple second-order 
models to represent appropriate underdamped and overdamped conditions. A similar 
procedure is followed to formulate the disturbance-rejection objectives of the design: In 
this case, the magnitude of the sensitivity function \S(p,jUi)\ = |(1 + G(p, jUi)K(joji))~ 1 \ 
is required not to exceed an appropriate bound D u (u{) for all % = 1,2,..., TV and every 
p G V. Note that the disturbance bounds are defined independently of the tracking 
bounds due to the presence of the prefilter, F(s). 

Next, the tracking and disturbance-rejection specifications are translated to certain 
conditions on the nominal open-loop frequency response L (ju>) = G (jcu)K(ju) where 
G (juj) denotes the nominal plant, defined for any p G V . Consider first the tracking 
bounds. For each frequency Ui, i = 1, 2, N, the tracking specifications will be satisfied 
if and only if, 

< 6(wi) := B u (ui)\ dh - Bi(cji)\ dh (1) 

db 

i.e. if and only if the maximum variation in the closed-loop gain as p varies over the set V, 
at each frequency Ui does not exceed the maximum allowable spread in the specifications 
8(uji). This is because, (i) the uncertainty associated with K is assumed negligible, and (ii) 
the actual closed- loop gain can always be adjusted to its required level at each frequency 
via the scaling action of the prefilter. For each frequency Wj, the open- loop gain (for 
each open- loop phase 0) at which condition (1) is met with equality, defines the so-called 
"Horowitz template" h(uji, (f>); this is the minimum open-loop gain necessary to achieve 
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max 



G(p,jUj)K(jUj) 
1 + G{p,jbJi)K(jui 



the specified robust tracking specifications at frequency out. In practice, each "Horowitz 
template" is calculated via a bisection algorithm over a phase grid, and the equality 
condition in (1) is satisfied within a given tolerance. In total, we have N Horowitz 
templates, one for each frequency of interest Ui. A similar procedure will result in N 
robust disturbance-rejection contours, di(ui, <fi), i = 1, 2, N, which define the minimum 
open-loop gain required to achieve robust disturbance rejection (i.e. disturbance rejection 
for every p G V). Clearly, to achieve robust tracking and robust disturbance-rejection 
simultaneously, the open-loop frequency response must satisfy, 

\L{jUi)\ > fi{ui,4>) := m&x{hi(ui,4>),di(ui,(f))} (2) 

for all i = 1, 2, N, where the maximum in (2) is calculated pointwise in = argL(jco>j). 
The contours fi(u)i, <p) will be referred to as the robust-performance bounds. 

In addition to robust-performance objectives, the closed-loop system should also be 
robustly stable, i.e. the closed-loop transfer function F(s)G(p, s)K(s)(l + G(p, s)K(s))~ 1 
should be stable for every p G V. Since F(s) will be designed stable, robust stability 
can be inferred from the number of encirclements around the —1 point by the open- 
loop frequency response L(p,ju) = G(p, ju)K(ju), p G V. In practice, a more severe 
constraint is imposed on L(p,ju>): To establish a minimum amount of damping for the 
(nominal) closed-loop system, the nominal open-loop response is constrained not to enter 
an M-circle of an appropriate value. Under the assumption of parametric uncertainty 1 in 
the plant, the uncertainty templates of G Q (ju) at high frequencies approach a vertical line 
in the Nichols chart. Hence, to ensure that the specified minimum amount of damping 
is maintained at high frequencies despite the presence of uncertainty, the M circle is 
translated downwards by a specific amount, resulting in the so-called "universal high 
frequency boundary" , which should not be penetrated by the nominal open-loop response 
[3|. This contour will be denoted by B in the sequel. 

2 Optimal Design of PID controllers 

After constructing the contours fi(u)i, <fi) and B on the Nichols' chart, QFT normally pro- 
ceeds with the design of the feedback controller K(s). This normally involves frequency- 
shaping of the nominal open- loop frequency response L {ju), so that it does not penetrate 
the the B contour and \L {juii)\ lie on or above the robust-performace bounds fi{joJi,4>) 
for each U{. At this stage, the designer normally follows a phase- lead/lag compensation 
approach which involves a considerable trial-and-error element, and can be cumbersome, 
especially if the number of specified frequencies is large and the specifications are tight. 
If the specifications can not be achieved, the design objectives are assumed to be infea- 
sible. In this case, the specifications are normally relaxed and the design is repeated. If 

1 The method can be modified to also take into account unstructured high frequency uncertainty. 
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the specifications are feasible, the best design is considered to be the one for which the 
specifications are met as tightly as possible. This is in order to avoid the possibility of 
"overdesigning" the system using unnecessarily large gains/bandwidth, which can result 
in measurement noise amplification and potential instability due to parasitics and high- 
frequency unmodelled dynamics [3] , jl] . A compromise between controller complexity and 
a "tight" design has to be made in many cases. 

In this section we present a simple algorithm for designing PID controllers which are 
optimal in the QFT sense. PID (or "three-term") controllers are widely used in industry 
because they are simple and easy to tune. Our algorithm can be used to provide an 
adequate QFT design, or as the first step for designing a more complex controller pQ. In 
addition, the algorithm can easily tackle a large number of constraints and can be also be 
applied to multivariable systems, using the standard QFT approach [7], [B]. 

A PID controller has a transfer function, 

k- 

K pid (s) = k p + -j + k d s (3) 

and is therefore completely defined by the three terms k p (proportional gain), hi (integral 
gain) and kd (derivative gain). Its frequency response is 

k- 

Kpid(juj) = k p -j—+ jk d u (4) 
tu 

Suppose that K pid (s) is used in cascade with an uncertain plant G(p, s). Then, the nominal 
open-loop system has frequency response L (ju) = G (joj)K pi(i (ju), where G a (s) denotes 
the nominal plant. Thus, the asymptotic gain of the nominal open loop system is given 
by, 



lim^oo 



G a (ju 




(5) 



Suppose that the asymptotic gain of the nominal plant is \G (ju))\ ~ Au~ p where the 
pole/zero excess of the nominal plant p is at least equal to 2. Then the asymptotic gain 
of the nominal open loop is \L (ju)\ ~ A\kd\uj~ p+1 . Since A and p are fixed, the nominal 
open- loop gain at high frequencies is minimised by minimising \ka\- This objective is 
consistent with the aims of QFT theory outlined previously. 

To design an optimal PID controller, consistent with the requirements of QFT, we 
formulate the following optimisation problem: 

Minimise \kj\ subject to the constraints: 

1. \L (jui)\ > fi(u h <p) for all % = 1, 2, N where <p := aigL (jui). 

2. L (jUi) £ B for all i = 1, 2, N. 

Since B is a closed contour which is defined for a range of phases only, it is always 
possible to combine constraints (1) and (2) for each i — 1, 2, ./V into a single constraint 
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as long as intersects B or lies entirely above it, by calculating the point-wise 

maximum of the two contours in the common phase range. This is almost always the 
case in practice, since robust performance objectives are almost never associated with 
frequencies significantly exceeding the closed-loop bandwidth. The (unlikely) case that a 
performance bound lies below the B contour can also be accomodated in our algorithm via 
an additional checking condition. To simplify the presentation, however, we will assume 
in the sequel that this does not occur, and the combined contours will be denoted by 
fi(u)i,(f>i). The optimisation problem, therefore, takes the form: Minimise \kd\ subject to 
Wu3i)\ > fi(u>i,<i>) for alH = l,2,...,iV. 

The magnitude (linear) and phase of K pid (ju) are given by, 



\K pid (ju)\ = 



k 2 v + (k d uj - ^ 



arg(K pid (ju)) := ip(u) = tan 1 



k n 



(6) 



respectively. Now suppose that we fix the phase of K pid at two distinct frequencies and 
u>j, i.e. ip(jcui) = ipi and ip(jujj) = ipj- Then, 

'k d Ui - £ 



which implies that 



Similarly, 



=ipi = tan 1 



ki /Cptan(^j) 



tot 



k d 2 



ki k p tan(ipj) 



kn 



= 



= 



(7) 



(8) 



(9) 



"3 ^3 

Equations (8) and (9) can be arranged in matrix form as: 



1 — 

t 

1 — 



tan(i/>j) 

OJi 

tan(^j) 



/ k d \ 

ki 

V K } 



(10) 



Denote the 2x3 matrix in equation (10) by A(t/j i: tpj). Clearly, Rank(A(-^j, ipj)) = 2 since 
OJi 7^ ojj. Therefore, the kernel of A(ipi,ipj) is a one-dimensional subspace of 7Z 3 , which 
implies that the controller gains kd, h and k p are fixed up to scaling. Numerically, the 
kernel of A(i/j i: i/jj) can be calculated easily using the Singular Value Decomposition. 
Applying the Singular Value Decomposition to A(t/j i: ipj) gives: 




(11) 



where V± is a 2 x 3 matrix. Here V 2 spans the kernel of A(ipi, ipj). Write = [V 2 i V 22 V23] 
Then, 



;i2) 



/ kd \ 




f V21 \ 


ki 


= A 




V K J 




{ V23 ) 
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where A is an arbitrary real constant. Using equation (6), the gain and phase of K pid (s) 
may be written as 

\K pid (ju)\ = 1X^,1+ (v 21 u-^\ = tan" 1 y~ ^ ) (13) 

Note that equation (13) implies that fixing the phase of K pid (ju) at two frequencies, fixes 
the phase of K p i d (juj) at any frequency u, and thus also the phase of L Q (ju). In this case, 
minimising \kd\ is equivalent to minimising | AV21 1 - 

Under the constraint that ipioJi) = ipi and ip(cjj) = ipj, the QFT constraints are 
satisfied if and only if 

\L Q (juj k )\ dh > f k (u k ,(f) k ) (14) 

for all k — 1,2, ...,N, where <p k = &rg(L (ju k )). Note that since the phase of K p \ d (joo) is 
fixed at every frequency u, <j) k is fixed and known for all k — 1, 2, N . In fact, 

, fV 2 itu k - ^\ 

<P k = exg(G (juj k )) + tan" 1 ^—^ (15) 

Since, 

\L {ju k )\ dh = \G {ju k )\ db + \K pid (ju k )\ dh (16) 
equation (14) is equivalent to 

|-^ P id(jwfc)|db > h(uk,<f>k) - \G (ju k )\ dh (17) 
for all k — 1,2, ...N. Substituting from (13) shows that this is equivalent to 



V22 



(18) 



201og 10 |A| > max fc= i i2 ,...,7v \ fk(wk,i/>k) - \G (ju k )\ dh - 101og 10 \ V 2 \ + (v 2 iu k 

or |A| > 10 20 where we have defined, 
= max fce{ i )2 ,...Ar} lf k (uj k ,(f) k ) - |G (ja; fe )|db - lOlogio \ V 23 + (v 21 u k - J J (19) 



Multiplying by | V21 1 and noting that \kd\ — | AV21 1 , implies that \k d \ > | V21 1 10 20 . Hence, 
provided that the phase of K pid is fixed as arg(i^ pi d(jco'i)) = ipi and aig{K pid {jojj)) = ipj, 
the minimum value of \k d \ which achieves the robust-performance constraints is given by 
\k* d \ = I V 2 fl0^~ I, where the additional indexes introduced in V 2 \ and (3 emphasise 
the dependence of these variables on (ui,Uj) and We can now formulate the 

following algorithm for solving the optimisation problem. In this algorithm, the gains 
ki, kd and k p have been further constrained to be non- negative. This assumption can be 
removed, if desired, with minor modifications to the algorithm. 

Algorithm 1: Given the plant's uncertainty templates Qi, a nominal plant G (jui) G 
Qi and QFT contraint bounds fi(uji,<f>), each defined at frequencies ui,u 2 , ...,un, the 
following algorithm calculates an optimal PID controller K* id (s) = k* + ^- + k* d s with 
non-negative gains, if it exists: 
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1. Obtain a phase array $ by discretising the phase interval (— 360°, 0°). 

2. Select any two distinct frequencies ui^Ui G {u>i, u 2 , ...,u>n}. 

3. Caclulate phase intervals <3>fc, C $ in which the nominal open-loop phase can vary 
at cu k and uj\ if a PID controller is used. $ fc ($;) contains the phases of $ which lie 
within ±90° of the nominal plant phase G (ju k ) (G (ui)). 

4. Initialise m x n arrays K p , Ki and K d , where m and n are the sizes of and $j 
respectively. 

5. For each ($ fc (i), G x 

(a) Calculate ^ = - arg(G (w fc )) and = $j0') - arg(G (<^)). 

(b) Calculate the singular value decomposition of A(ipi, ip k ) and the corresponding 
vector V 2 3 spanning its kernel. 

(c) If any two elements of V 2 3 have opposite signs, set K d (i,j) = 00; Else, let q be 
the (common) sign of the elements of V^ 1 , calculate /3 i: j using equation (17), 
and set A = qlO^ , K d (i,j) = \V 2 \ j , K&J) = \V 2 % j and K p (i,j) = \V 2 \ j . 

6. Calculate € aigmin i: jKd(i, j). If K d {i*,j*) = 00, the constraints cannot be 
satisfied by a PID controller with non-negative gains; otherwise the optimal PID 
controller is K*(s) = K p (i*,f) + K^i*, j*)s _1 + K d (i*,j*)s. 

Remarks on Algorithm 1: 

• In step 1 the phase discretisation of the interval (— 360°, 0°) results in a grid of 
phases, typically equally-spaced. In practice, 50-100 phases are adequate. It is 
helpful to calculate the performance bounds over the same phase grid (<£>). 

• In principle any two frequencies u k , io\ can be selected in step 2. Selecting these two 
frequencies reasonably far apart works well in practice. 

• In step 3, the phase intervals and $; may be further restricted, if desired, to en- 
sure that the nominal open-loop frequency response L Q (ju) is shaped appropriately. 
This will also reduce the number of calculations in step 5 of the algorithm. 

• The singular value decomposition in step 5(b) of the algorithm can be dispensed with 
altogether, by calculating V 2 1 ' 3 analytically. This, however, increases the complexity 
of the algorithm and does not lead to any significant reduction in computation time. 

• Step 5(c) of the algorithm requires the calculation of flij which in turn relies on the 
calculation of the performance bounds at phases specified by equation (15), which 
may not belong to the (discrete) phase array used to calculate the bounds (typically 
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$). There is no difficulty, however, in estimating the performance bounds for these 
phases using the adjacent phase points, e.g. via linear interpolation. Alternatively, 
the performance bounds can be calculated exactly for the phases obtained from 
equation (15) via a bisection algorithm implemented between steps 5(a) and 5(c). 

• The algorithm can be specialised, if desired, to calculate optimal PI of PD con- 
trollers. This requires only a one-dimensional search over a phase grid defined at a 
single frequency. 



3 Example 

In this section we illustrate our algorithm by means of a simple example. The uncertain 
plant is taken as 

G(s) = J^- (20) 
s z + as 

in which the parameters a and k vary independently in the intervals 1 < a < 10 and 
1 < k < 10 respectively. The nominal plant, G Q (s), is taken to correspond to a = 1 and 
k = 1. The tracking specifications are defined as: 

\Bi(jUi)\ < \T(a,k,jUi)\ < \B u (jui)\ (21) 

Here T(a, k, s) = F(s)G(a, k, s)K(s)(l+G(a, k, s)K(s))~ 1 is the closed-loop transfer func- 
tion and the lower and upper tracking bounds are defined as the magnitude frequency 
responses of the two systems 

B,( S ) = 7 °- 658 t + 30) r (22) 

KJ (s + 2 + j3.969)(s + 2- j'3.969) v ; 

and 8400 

Bu{s) = ( s + 3)( S + 4)( S + 10)( S + 70) (23) 
at each s = ju{. Note that the zero of Bi(s) (s = —30) and the two fast poles in B u (s) 
(s = —10 and s = —70) have been included to ensure that the the magnitude frequency 
responses of Bi(s) and B u (s) diverge at high frequencies [2]. The frequencies of interest 
{oji} have been selected as u>i = 0.5, u 2 = 1, u 3 = 2, co> 4 = 3, co 5 = 5, ujq = 10, u 7 = 30 and 
us = 60 rads/s. For simplicity, no disturbance-rejection objectives have been considered 
in this example. 

The uncertainty templates of the plant at the eight frequencies of interest have been 
calculated numerically and are displayed in Figure 2. To reduce the number of calcu- 
lations, each template has been replaced by its convex hull. This results in a minimal 
amount of conservativeness in this case. 
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Figure 2: Uncertainty templates 

Next, an optimal PID controller was designed following the procedure of Algorithm 1. 
The optimal controller was obtained as 

4.46 

K; id (s) = 12.6 + 3.95s + (24) 

s 

Figure 3 shows the frequency response of the nominal plant (dashed line) and the 
nominal open loop (solid line) on the Nichols chart, together with the eight Horowitz 
templates and the £>-contour (corresponding to an M value of 1.2). The eight frequencies 
of interest are indicated by circles on the two frequency responses. The design meets the 
specifications, since the nominal frequency response does not penetrate the S-contour and 
lies on or above the Horowitz templates at the eight frequencies of interest. As expected, 
one of the L(jo;j)'s (L(ju^)) lies exactly on a bound (in this case the B contour). 

Since the open-loop system has a pole-zero excess equal to one, its phase approaches 
—90° at high frequencies. There is no difficulty, however, in forcing the open loop response 
to approach the —180° phase line at high frequencies, if desired, by including a suitable 
1 + sr factor in the denominator of the controller derivative term. Choosing, for example, 
r <C Wjy 1 , has a minimal effect on the shape of L (juu) for uj < ujn- Alternatively, the PID 
controller may be assumed to be of the form 

= + 7 + IT^ (25) 
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Figure 3: Nominal plant and nominal open-loop 

(r fixed) before solving the optimisation problem. Since in this case 

L' (p, s) := K' pid (s)G(p, s) = ((k p + hr) + ^ + (k d + k p r)s) (26) 

y s J 1 + st 

our algorithm can still be applied by redefining the uncertain plant as 

0-fe.) = ^f> (27) 

1 + ST 

and optimising with respect to the new variables k' = k p + k{T, k[ = ki and k' d = kd + k p T. 

A second-order prefilter F(s) (of dc gain equal to 1 and cut-off frequencies 3.5 and 
7.5 rads/s) was finally designed using the standard procedure [2]. Figure 4 shows the 
closed-loop frequency responses for a number of (a, k) parameter combinations, together 
with the tracking bounds \Bi(joj)\ and \B u (jcu)\. Again, the eight frequencies of interest 
are marked by circles. Clearly, the specifications of the design are met, as expected from 
the characteristics of the open loop response in Figure 3. 



4 Conclusions 

An algorithm has been presented for designing optimal PID controllers for uncertain 
systems subject to QFT constraints. The algorithm is simple, easy to implement and can 
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Figure 4: Closed-loop responses and bounds 

be used to automate the loop-shaping step of the QFT design procedure. Although the 
algorithm has been presented for SISO systems, its extension to multivariable problems 
is possible using the standard QFT approach. 
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